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In a previous paper [gr-qc/0104001 ; Class. Quant. Grav. 18 (2001) 3595-3610] we have shown 
that the occurrence of curved spacetime "effective Lorentzian geometries" is a generic result of 
linearizing an arbitrary classical field theory around some non-trivial background configuration. 
This observation explains the ubiquitous nature of the "analog models" for general relativity that 
have recently been developed based on condensed matter physics. In the simple (single scalar field) 
situation analyzed in our previous paper, there is a single unique effective metric; more complicated 
situations can lead to bi-metric and multi-metric theories. In the present paper we will investigate 
the conditions required to keep the situation under control and compatible with experiment — 
either by enforcing a unique effective metric (as would be required to be strictly compatible with 
the Einstein Equivalence Principle), or at the worst by arranging things so that there are multiple 
metrics that are all "close" to each other (in order to be compatible with the Eotvos experiment). 
The algebraically most general situation leads to a physical model whose mathematical description 
requires an extension of the usual notion of Finsler geometry to a Lorentzian-signature pseudo- 
Finsler geometry; while this is possibly of some interest in its own right, this particular case does 
not seem to be immediately relevant for either particle physics or gravitation. The key result is that 
wide classes of theories lend themselves to an effective metric description. This observation provides 
further evidence that the notion of "analog gravity" is rather generic. 
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I. INTRODUCTION 



Whenever you have a single classical scalar field governed by Lagrangian dynamics, the linearization of the field 
around any classical solution automatically provides you with a curved space "effective metric" . Furthermore this 
effective metric has Lorentzian signature if and only if the PDE governing the scalar field is hyperbolic This 
observation, combined with the fact that quantization of the linearized fluctuations automatically provides an Einstein- 
Hilbert term in the one-loop effective action, led us to suggest that Einstein gravity is a generic low-energy limit of a 
wide class of quantum field theories [ij. 

A technical step in |0] that kept the calculation under control was to assume that one was dealing with a single 
scalar field, which automatically implied that there was a single unique effective metric. However, this assumption 
is very restrictive. By starting with a single scalar field it is not possible to reproduce the entire algebraic structure 
of the set of possible metrics in GR. Therefore, in the present article we turn our attention towards the somewhat 
messier general situation of multiple fields (multiple scalar fields, or a multi-component vector or tensor field). The 
key new results in this situation are: 

1. In some situations the existence of a single unique effective metric can be deduced. This is the case that is 
compatible with strict application of the Einstein Equivalence Principle. 
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2. More generally there may be several distinct effective metrics on the same spacetime (this will generically occur 
if the dynamical system exhibits "birefringence" — more properly called "multi-refringence" — a situation in 
which different types of oscillation, typically referred to as different "polarizations" , can propagate at different 
speed). As long as these various metrics are sufficiently "close" to each other, different types of matter will feel 
approximately the same geometry — as is required experimentally to satisfy the constraints deduced from the 
Eotvos experiment. 

3. In the algebraically most general situation a mathematical structure is encountered which can best be viewed 
as a hyperbolic extension of the notion of a Finsler geometry; these pseudo-Finsler geometries bear the same 
relation to ordinary Finsler geometries that pseudo-Riemannian geometries {aka Lorentzian geometries) bear 
to Riemannian geometries — and many of the same sort of technical problems arise due to indefiniteness of the 
pseudo-Finslerian "metric". While these structures may be of interest in their own right, they do not seem to 
be immediately appropriate for either particle physics or gravitation. 

In our previous paper we studied the single-scalar case in detail; that case being a perfect exemplar for the occurrence 
of a Lorentzian effective metric . That paper also contained an extensive bibliography regarding analog models and 
we will be more selective in this present paper. We shall consider multi-component systems and push the "effective 
metric" analysis as far as practical, making extensive use of the theory of characteristic surfaces. We derive algebraic 
conditions that should be satisfied to keep the effective metric unique (or nearly so). We investigate the notion of 
"polarization" and the associated "Fresnel equation" , leading to the notion of "multi-refringence" . Finally we have a 
few words to say about the algebraically general pseudo-Finsler geometries — this is a subject that mathematically is 
very poorly developed, and we hope by this paper to generate some interest in what otherwise seems a rather arcane 
and abstract subject. The key message to take from the present effort is that while dealing with multiple fields is 
algebraically more tricky, there are nevertheless wide classes of dynamical systems that lend themselves to an effective 
metric description. Furthermore with multiple background fields, the background geometry is more flexible — the 
metric could in principle be completely general. 



II. LINEARIZED FIELDS FROM ARBITRARY BACKGROUND SYSTEMS 



A. Lagrangian analysis 



Suppose we consider a collection of fields — {(f>^ , (f)^ , . . .} whose dynamics is governed by some first-order 



Lagrangian C{df_i 



). Here "first-order" is taken to mean that the Lagrangian depends only on the fields and 



their first derivatives. We want to consider linearized fluctuations around some background solution of the equations 
of motion. As in the single-field case [|| we write 



(1) 



Now use this to expand the Lagrangian 
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Consider the action 



+ 0{e'). (2) 



(3) 



Doing so allows us to integrate by parts. As in the single-field case we can use the Euler-Lagrange equations to 
discard the linear terms (since we are linearizing around a solution of the equations of motion) and so get 
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Because the fields now carry indices (AB) we cannot cast the action into quite as simple a form as was possible in 
the single-field case. The equation of motion for the linearized fluctuations are now read off as 



This is a linear second-order system of partial differential equations with position-dependent coefficients. This system 
of PDEs is automatically self-adjoint (with respect to the trivial "flat" measure d'^+^x). 
To simplify the notation we introduce a number of definitions. First 

This quantity is independently symmetric under interchange of ^, v and A, B. In the case of Bose fields we will want 
to interpret this as something related to some sort of "metric" , but the interpretation is not as straightforward as for 
the single-field case. 
Next, define 

^ d^C d^C 1 / d^C d^C \ 

^9(9^0^) did^(j)^) d^"^ ^ 2 " \did,<j)'^) did^(j)^) 9(9^0^) 9(5.0^)7 ■ ^' 

This quantity is anti-symmetric in A, B. For Bose fields we will want to interpret this as some sort of "connexion". 
Equivalently we could write 

^9(9^0-4) 9<^s 9(9^<^S)5^A 2 \d{d^<t^^) d{d,<t)B) d{d^<l>B) d{d,cj)^) 



It is useful to note that 

In the case of Fermi fields one typically has J^'^ab = 0, while the T^ab arc most usefully thought of as generalizations 
of the Dirac matrices. Finally, define 

d'^C 1 / 9^£ \ ^ f \ 
^^"^ " "90^ d^B + \d{d^^^) 90S j ^ 2^^ \d{d^cj>B) 90-4 j ■ (10) 

This quantity is by construction symmetric in {AB). We will want to interpret this as some sort of "potential" or 
"mass matrix" . Then the crucial point for the following discussion is to realize that equation (|^) can be written in 
the form 

9^ {rAB 9.0f ) + \ [T% d^cj^f + 9,,(r^5 0f )] + Kab cj>f = 0. (11) 
Now it is more transparent that this is a formally self-adjoint second-order linear system of PDEs. 



B. Systems of second-order hyperbolic PDEs 

One can arrive at a similar set of equations by directly linearizing an arbitrary system of second-order PDEs, that is, 
without relying on the existence of a Lagrangian dynamics. Start by considering an arbitrary system of second-order 
PDEs written in the form 

Ga{x, g^g^^S) ^ 0_ (^2) 

The PDE does not have to be linear or even quasi-linear; there are as many equations as there are fields. As in the 
single-field case defining hyperbolicity for such a general equation is not easy — not even Courant and Hilbert ||] 
consider this case explicitly. Indeed we shall adapt their discussion and shall define hyperbolicity in terms of the 
linearized equation: Suppose we linearize around some solution 0;^, writing 



(13) 
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Then 



dG, 



0. 



(14) 



That is: the fluctuation satisfies a second-order system of linear PDEs with time-dependent and position-dependent 
coefficients (these coefficients depend on the background field you are linearizing around). This equation can be 
written in a manner similar to ( [ll|) 



by defining 



(r^AB 9.0f + d^(t>f + Kab 0f = 0, 
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(16) 
(17) 
(18) 



Though the form of the coefRcients is somewhat different from that appearing in the Lagrangian-based analysis (in 
general the equation is not automatically self-adjoint), much of the discussion that follows can be applied (with minor 
modifications) to both cases. 



C. Eikonal approximation 



The theory of systems of second-order PDE is relatively complicated and much less transparent than that for a 
single second-order PDE. (See, for example, Courant and Hilbert j^, volume 2, pp 577-618, or the Encyclopedic 
Dictionary of Mathematics |^.) We would like to be able to construct some notion of spacetime metric directly from 
the quantities f^'^AB (or J^'^ab), but we shall soon see that doing so is a somewhat tricky proposition. 

Consider an eikonal approximation for an arbitrary direction in field space, that is, take 

0^(a;) = e^{x) exp[-*^(a;)], (19) 

with e^{x) a slowly varying amplitude, and '■p{x) a rapidly varying phase. In this eikonal approximation (where we 
neglect gradients in the amplitude, and gradients in the coefRcients of the PDEs, retaining only the gradients of the 
phase) the linearized system of PDEs ( |ll| ) becomes 

{f'^AB d^^^ix) dM^) + ^^AB dM^) + Kab} ef = 0. (20) 

This has a nontrivial solution if and only if e^{x) is a null eigenvector of the matrix 

f'^'AB d^x) dM^) + r^AB df^^ix) + Kab- (21) 

Now, the condition for such a null eigenvector to exist is that 

Fip, q) EE det {r^AB dM^) dM^) + r^AB dM^) + Kab} = 0, (22) 

with the determinant to be taken on the field space indices. This is the natural generalization to the current situation 
of the Fresnel equation of bi-refringent optics ^ ^ . Following the analogy with the situation in electrodynamics (either 
nonlinear electrodynamics, or more prosaically propagation in a bi-refringent crystal), the null eigenvector €^(x) would 
correspond to a specific "polarization" . The Fresnel equation then describes how different polarizations can propagate 
at different velocities (or in the language to be used later in the paper, can see different metric structures). In particle 
physics language this determinant condition F{p,q) = is the natural generalization of the "mass shell" constraint. 
Indeed it is useful to define the mass shell as a subset of the cotangent space by 

F{p,q) = o}. (23) 
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In more mathematical language we are looking at the null space of the determinant of the "symbol" of the system 
of PDEs. By investigating F{p, q) one can recover part (not all) of the information encoded in the matrices J^'^ab, 
Tf^AB, and Kab, or equivalently in the "generalized Fresnel equation" (|2^). (Note that for the determinant equation 
to be useful it should be non- vacuous; in particular one should carefully eliminate all gauge and spurious degrees of 
freedom before constructing this "generalized Fresnel equation" , since otherwise the deter minan t will be identically 



zero. Some examples of this phenomenon will be given later in the paper, see subsection IVE ). We now want to 
make this analogy with optics more precise, by carefully considering the notion of characteristics and characteristic 
surfaces. We will see how to extract from the the high-frequency high-momentum regime described by the eikonal 
approximation all the information concerning the causal structure of the theory. 

III. CAUSAL STRUCTURES 

One of the key structures that a Lorentzian spacetime metric provides is the notion of causal relationships. This 
suggests that it may be profitable to try to work backwards from the causal structure to determine a Lorentzian metric. 
Now the causal structure implicit in the system of second-order PDEs given in equation ( |ll|) is described in terms of 
the characteristic surfaces, and it is for this reason that we now focus on characteristics as a way of encoding causal 
structure, and as a surrogate for some notion of Lorentzian metric. Note that via the Hadamard theory of surfaces of 
discontinuity the characteristics can be identified with the infinite-momentum limit of the eikonal approximation 
That is, when extracting the characteristic surfaces we neglect subdominant terms in the generalized Fresnel equation 
and focus only on the leading term in the symbol {J^'^ab)- In particle physics language going to the infinite-momentum 
limit puts us on the light cone instead of the mass shell; and it is the light cone that is more useful in determining 
causal structure. The "normal cone" at some specified point q, consisting of the locus of normals to the characteristic 
surfaces, is defined by 

A/'(9) = (pm 



det (/'"'^s p^) - 0| . (24) 

As was the case for the Fresnel equation (^2|), the determinant is to be taken on the field indices AB. Remember to 
eliminate spurious and gauge degrees of freedom so that this determinant is not identically zero. (See Courant and 
Hilbert ||^, volume 2 "Partial differential equations" , page 580.) We emphasise that the algebraic equation defining the 
normal cone is the leading term in the Fresnel equation encountered in discussing the eikonal approximation. If there 
are N fields in total then this "normal cone" will generically consist of nested sheets each with the topology (not 
necessarily the geometry) of a cone. Often several of these cones will coincide, which is not particularly troublesome, 
but unfortunately it is also common for some of these cones to be degenerate, which is more problematic. See Courant 
and Hilbert footnote 1 on page 592: 

It may he remarked that the present state of the theory of algebraic surfaces does not permit entirely 
satisfactory applications to the questions of reality of geometric structures which confront us here. 

Note that if one is dealing with Fermi fields where J^'^ab = 0, one has a choice: either iterate the first-order system 
of PDEs to produce a second-order system and then apply the previous logic, or go to the highest nontrivial term in 
the symbol and rephrase the discussion below [mutatis mutandis) in terms of 



det [T^ab Pm) = 0| • (25) 
Returning to the Bosonic case, it is convenient to define a function Q{q,p) on the co-tangent bundle 

Qiq,p) = det {r''AB{q)p^P^.). (26) 

The function Q{q,p) defines a completely-symmetric spacetime tensor (actually, a tensor density) with 2N indices 

Q(g,p) = QMi-iM.-2-p«.«(^) p^^... p^^ p^^_ (27) 

(Remember that J^'^ ab is symmetric in both and AB independently.) Explicitly, using the expansion of the 
determinant in terms of completely antisymmetric field-space Levi-Civita tensors 

Q^lVl^2V2 - tJ.NVN _ ^AiA2A3---An ^BiB2B3---Bn (28) 

' F'^^'mB, f'-'A^B.^ 
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In terms of this Q{q,p) function the normal cone is 



Q{q,p)^0\. (29) 



In contrast, the "Monge cone" {aka "ray cone", aka "characteristic cone", aka "null cone") is the envelope of the 
set of characteristic surfaces through the point q. Thus the "Monge cone" is dual to the "normal cone" , its explicit 
construction is given by (Courant and Hilbert ||], volume 2, pp 583): 



Miq)^{t^^^^ 
op,. 



p,, e U{q) . (30) 



Equivalently 

M{q) = H'^ = Q'"^--''-(q) p^, p^, ■■■p^. 



e ma) ■ (31) 



Unfortunately, even if the normal cone is pleasantly behaved, the Monge cone (null cone) may be quite messy. Indeed 
Courant and Hilbert remark: 

Even if [the normal cone] is a relatively simple algebraic cone of degree [2N], the ray cone [Monge cone/null 
cone] may have singularities, or isolated rays, and need not consist of separate smooth conical shells. 

(Courant and Hilbert volume 2, p 584). Defining the notion of "hyperbolicity" in this general context boils down 
to the question of just what constraints have to be placed on the original system of PDEs to make sure the Monge 
cone is well-behaved. 

For completeness, we mention a general construction for finding bi-characteristic rays (these are lines which are 
guaranteed to lie on characteristic surfaces). Start by picking some point q and a co-tangent vector p which is in the 
normal cone M{q). Now with these initial conditions solve the "Hamiltonian equation" 



dqf" ^ dQ{q,p) ^ dp^ ^ ^ dQ{q,p) 
ds ' ds dq>^ 



(32) 



Here s is just a parameter, it is not physical time. The resulting curve [q'^{s), p,i{s)] is called a bi-characteristic ray. 
For all s we have p,i{s) G Af{q{s)) and dq'^/ds E A4{q{s)). (The momentum lies on the normal cone, and the velocity 
lies on the Monge cone [null cone].) 

In a manner similar to the way in which the Monge cone is the dual of the normal cone, we can define a dual to 
the mass shell, this dual now being a subset of the tangent space 



T*{q) 



f dF{p,q) 
I dp,. 



P e Hq) ■ (33) 



In simple situations where there is a single metric the Monge cone can be constructed simply by "raising the index" 
in the definition of the normal cone, and similarly for the dual mass shell in terms of the mass shell. 

The structure of the normal and Monge cones encode all the information related with the causal propagation of 
signals associated with the system of PDEs. In next section, we will detail how to relate this causal structure with 
the existence of some spacetime metrics in different specific situations, from the experimentally favored single-metric 
theory deduced from the Equivalence Principle to the most complicated case of pseudo-Finsler geometries. 

IV. GEOMETRICAL INTERPRETATION 
A. Field redefinitions 

We have seen that the causal structure of a system of coupled PDEs of the form (^ could in general be rather 
complicated. However, there are particular cases in which it is relatively easy to find a geometrical interpretation of 
what is happening. First, it is important to realize that we are always free to perform a field redefinition 

0^ ^ 0^ = /i^(0S), (34) 
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essentially a coordinate change in field space. This induces a linear transformation on the linearized fields (p^, 



dh 



A 



f (35) 



and, therefore, a redefinition of the f^'^ ab, ^^cd, and Kab- It is convenient to adopt matrix notation, suppressing 
the field indices (but keeping the spacetime indices explicit). In matrix notation 

01 ^ 01 - i-i 01, (36) 

while the self-adjoint system of PDEs ( pi] ) can be written in the form 

1 



d,cl>i) + 9^01 + -d^{T^) (Pi+Kcf,i^O. (37) 



A brief computation yields 

f"" f"" = f"" L; (38) 

^ f ^ = i + f"" d^L - d^L^ f"" L; (39) 

K ^ K = L^ K L ]^ {duL^ V L V d^L) + {L^ f"" d^L + d^L^ /^"^ L) 

-d^L^ d^L. (40) 

Note that because of the matrix multiplications the ordering is important. Also note that {f^'^Y' — /^''i while 
(r'^)^ = — r'' and {Ky^ = K, where the transpose is on the field indices (not the spacetime indices), and that 
these symmetry properties are preserved under these background-field-dependent coordinate transformations. The 
point is that we can use these field redefinitions to simplify the matrix J^'^ab using position-dependent general linear 
transformations — this will be useful below. (For our purposes, the detailed form of the non-homogeneous derivative 
terms is often not important, though the fact of their presence is.) 

B. Einstein Equivalence Principle 

The physically simplest situation one can encounter is when the multiple fields being analyzed simply co rresp ond 



one 



to different but equivalent polarizations of a single field. By looking at the eikonal approximation of section 11 A, 
can verify if there exist multiple independent polarizations that are solutions of equation ( pO| ) for the same phase 
function Lp. Those multiple solutions should be grouped into a single class, all having the same propagation features. 
The number of independent solutions grouped into equivalent classes will give us information about the degree of 
degeneracy of the matrix /^"ab- (Here we arc thinking of this object as a matrix in field space, whose elements are 
spacetime matrices). 

The simplest possible situation regarding the geometrical structure emerging from a field theoretical mode analysis 
is that in which all linearized fields can be grouped into a single class. This implies that there must be some choice 
of field variables so that all the new (jyf see the same metric, that is: 

rAB = 5ab Z^" - 6ab g^r (41) 

This corresponds to the field theoretical analog PDE systems obeying strict adherence to the Einstein Equivalence 
Principle. If we use any other choice of field variables then we must have 

f'^AB = hAB = hAB (42) 

This "factorization" condition on f^'^AB is a necessary and sufficient condition for strict adherence to the Einstein 
Equivalence Principle. Indeed for strict adherence to the Einstein Equivalence Principle you would also want to 
demand that in the 0^ field variables ab — 0, and one would want the "mass matrix" Kab to be position 
independent. (At the very worst Kab might contain curvature coupling terms that would go to zero in the fiat space 
limit.) 

In the usual formulation of general relativity, strict adherence to the Einstein Equivalence Principle is enforced by 
a policy of "minimal substitution" — the matter Lagrangian is taken to be a flat Minkowski-space Lagrangian with 
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the substitution ry^j, g^v Since, by fiat, there is only one metric put into this Lagrangian, a unique spacetime 
metric will emerge for all the matter fields. (At least as long as all "external fields" are set to zero.) While the 
Einstein Equivalence Principle is certainly compatible with experiment, strictly enforcing it may be overkill — indeed 
the low-energy field theories arising from string theory and its relatives (in particular, various scalar tensor theories) 
do not strictly adhere to the Einstein Equivalence Principle, so it is useful to keep a little flexibility on this point. 

C. Multiple metrics 

Suppose now that we are looking for a multi-metric theory, then there must be some choice of field variables so that 
all the linearized fields "decouple" and see independent metrics. That is, there must be a "diagonal" representation 
in field space such that 

= diag{/r^ /r, /3^^ • • • m = diag{v^ ^r, 52^^ ■ ■ ■ .9^}- (43) 

If we use any other choice of field variables then the necessary and sufficient condition for the f^'^ ab to be simulta- 
neously diagonalizable in field space is that V /z, v, a, (3: 

f^AB r^Bc^r^AB f'"'BC- that is [ f^" , ] = 0. (44) 

Subject to this constraint, after diagonalizing f^'^AB in field space, each field couples to its own independent spacetime 
metric. If the original system of PDEs [equation (|ll|)] is hyperbolic then all these metrics will have Lorentzian 
signature. Note that this is the behaviour that typically occurs in nonlinear electrodynamics and in birefringent 
crystals; and note that such birefringence is not in conflict with the physical intent of the Einstein Equivalence 
Principle because in both of these situations there is an external field (the electromagnetic field; the rest frame 
of the crystal itself) that softly breaks the Lorentz invariance. One thing that we can definitely say based on the 
Eotvos experiment [the observed universality of free fall] is that those low-energy quantum fields relevant to describing 
ordinary bulk matter all see approximately the same effective metric. 

Note that the case of multiple metrics as defined above is related to (though not identical to) the notion of 
"reducibility" as introduced by Courant and Hilbert p 596. The function Q{q,p) is said to be "reducible" when 
it factorizes into lower-order polynomials. If it factorizes completely into simple quadratic products then: 

N N 

Q{q,p) = n QA{q,P) = n ^^a{<i) P^^ ■ (45) 

A=l A=l 

If this "quadratic reducibility" property is satisfied then the normal cone consists of N nested topological cones each 
of which is geometrically a cone, and we can define a variant "multiple-metric theory" by defining the N spacetime 
metrics using the prescription 

97 = (46) 

Even then, in order for the system ( pT| ) to be deemed to be hyperbolic, one must additionally insist that each of 
these matrices have Lorentzian signature and that the causal structures derived from each of these Lorentzian metric 
be compatible with each other: not only should each Lorentzian metric satisfy some sort of "chronology condition" 
or "causality condition" (no closed null or timelike paths), but to prevent an ill-posed problem one must insist that 
one cannot even form closed causal loops by using a chain of causal line segments belonging to different Lorentzian 
metrics. 

Now if we have a multiple metric theory defined in terms of simultaneously diagonalizing the ab this automat- 
ically satisfies quadratic reducibility for Q(q,p). The converse is not necessarily true, and the class of multiple-metric 
theories defined in terms of quadratic reducibility is more general than that defined in terms of simultaneous diago- 
nalization of the kinetic terms. (This last comment is not supposed to be obvious; we shall give an explicit example 
of this phenomenon when we discuss the generic two-field situation in section 0.) 

In many cases all these metrics will coincide, and we can speak of the spacetime metric; this happens for instance 
(modulo some technical issues such as gauge fixing) for both the Maxwell equation and the Dirac equation. In terms 
oi Q{q,p) this requires 

Q{q,p) - [Qa{q,p)f = [^(9) P^f ■ (47) 
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D. The general case: pseudo-Finsler geometries 



The most general case, when J^'^ab cannot be described in terms of a single metric, or even multiple metrics, must 
be dealt with by the formalism of pseudo-Finsler geometries. What we call pseudo-Finsler geometries bear the same 
relation to Finsler geometries that pseudo-Riemannian geometries {aka Lorentzian geometries) bear to Riemannian 
geometries — and they exhibit many of the same sort of technical problems that arise due to indefiniteness of the 
"metric" . 

Remember that for Riemannian geometries all the physically interesting quantities (metric, Riemann tensor, etc.) 
can be encoded in terms of coincidence limits of the (real) distance function d{x,y) and its derivatives. If we try to 
extend this sort of analysis to pseudo-Riemannian geometries then using the distance function is awkward because it 
is sometimes real (spacelike separation), sometimes zero (null separation), and sometimes pure imaginary (timelike 
separation). Synge (and many others) have argued that the best way of taking care of this technical difficulty is to 
define the "world function" Q 

n{x,y) = ^d{x,yf (48) 

which is always guaranteed to be real (though it can be positive, zero, or negative). 

In the present situation we have found it useful to use the notion of characteristic surfaces to encode the causal 
structure of the system of PDEs given in equation (^) and to define the real quantity Q{q,p) which can be positive, 
zero, or negative. Now Q{q,p) is not homogeneous linear in p, rather it is homogeneous of order 2N: 

Q{q,Xp) = X^^Q{q,p). (49) 

If we try to define the analog of the normal Finsler distance function, we would choose 

dF{q,p) = [Q{q,p)]'/"'. (50) 

If Q{q,p) were always positive (which would correspond to an elliptic system of PDEs) this construction can be 
used to provide a real and positive distance function suitable for defining a Finsler geometry. Unfortunately, we 
are interested in the hyperbolic case, so Q{q,p) by definition possesses zeros and changes sign. So dpiq^p) is now 
generically complex; one typically encounters various branch cuts involving 2N'th roots of —1, (exp{z7r/(2A^)}), which 
invalidates the standard presentation of Finsler geometries. (We emphasize that this is a technical issue, not a 
fundamental issue, but it does mean one cannot simply copy results from the standard mathematics literature.) Note 
that even in the case of a unique spacetime metric (completely reducible, so that Q{q,p) = [Qo{q,p)]'^) one still 
has dp{q,p) = y/Qo{q,p) = do{q,p). In this case the Finsler function degenerates to the usual Lorentzian distance 
distance function (which is positive real, zero, or pure imaginary). In view of the above, we see that it is Q{p,q) 
itself that should be thought of as fundamental: it is the natural pseudo-Finslerian generalization of Syngc's world 
function. 

There have been several attempts at defining Lorentzian-signature pseudo-Finsler geometries (for example, 
Asanov jsj). Unfortunately those pseudo-Finsler geometries are typically set up in such a way as to avoid the 
possibility of multiple light cones, which is exactly the situation we are trying to probe in this article. It does not 
seem to us that the Asanov formulation of pseudo-Finsler geometries has anything to say about the situation at hand. 
While it is clear that in the general case we will want to invoke some form of pseudo-Finsler geometry, none of the 
extant formalisms are really suitable |lO[ ^ |l^ . We hope that these brief comments might stimulate some 
interest in further developing this field. 



E. Hidden geometries — Eliminating spurious fields 

There are cases in which straightforward application of the previous analysis of characteristics does not work 
because the generalized Fresnel equation ( p2|) is identically zero. This happens in situations in which one of the fields 
is spurious, either because of a gauge invariance or possibly because of some non-obvious algebraic (not differential) 
relation between the fields. However, it should be noted that in some of these cases one can eliminate the spurious 
field completely, and thereby show that the remaining physical fields are coupled to a "reduced" effective metric. 

For example, let us take the general set of equations (pi]). Consider the situation in which for a particular equation, 
say equation number 1, and a particular field, say 0}, we have f^^n — 0, and F'^n = 0, but Ku ^ 0. Now choose 
the notation a, b to denote those indices A, B distinct from 1. Then, the relevant first equation from the system ( pT| ) 
reads: 

Ku 01 + (/(T + i (r^, d^cf^l + d,, [v;, 4]) + k,, <j,\ = o. (51) 
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This means that we can algebraically solve for the field using 



(52) 



This can now be used to eliminate from the system ( |ll| ) leading to a reduced [{N — 1) x [N — 1)] system of 
equations. To do this we first write the remaining a equations as 



(/aT 5.0?) (/rr d^^ct^l) + \ {Kb + d,K, 4>\]) + \ + d,K, 4>\])+K,b 4+Kai cj^l = 0, (53) 

and then substitute 

In general, this reduced set of equations has up to fourth-order derivatives of the linearized fields and so cannot 
be analyzed along the lines developed in this paper. However, there exist some particular situations in which the 
reduced set of equations is still second order, though no longer formally self-adjoint. Specifically, within this method 
of eliminating spurious degrees of freedom, there are three possible particular cases: 



lb = and r^',, = 0, 



• Case a: f^"'^ 

• Case b: f'^'^ai = and T^^ = 0, 



• Case c: f'^u = and /' 



al 



0. 



In all three of these cases the reduced set of equations can be written in second-order form 



ab 



Kab = 0. 



(54) 



Specifically, we have: 
• Case a: 



Jab — J J al 



Kb = Kb-r''aid, 



K, 



ab 



Kib 


)-Ki 


Kib 












Ku 


Kii 














K,b\ 
KiJ_ 


-Kid, 


fKu\ 
KKiiJ 


~lid,Ki) 


Ku 
Kii " 


^^(^.r::,) 


- Kal -TT- 
J^ll 



(55) 
(56) 
(57) 



Case b: 



Jab 


J a 


fuu Kal 
b — J lb -T7 — , 
i^ll 


^ ab 


= Kb^ 




Kab 


= Kab 





(58) 
(59) 
(60) 



• Case c: 



Jab 



^ ab 



ah 



2K, 



(KiKb + KiKb)^ 



^ ab 



^i9.Ki)Kb 



2K, 



-Kiid.Kb) + d, 



2K, 



KiKb- 



2K, 



al^ lb 



Kab = 



Ku 
Kn 

Kab- 

Ku 
2K,. 



Kal 

Kn 



J- lb 



-T^ d 

2 ab^f^ 



J^ll 

K 



Kid, 



K 



^)-77^(«i)(W) 



{d,Ki) ~ ^(^mF^'J - l^KuT,{d,Kb)- 



Kn 

Kal 1 
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(61) 



(62) 



(63) 
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If one now forgets about the initial non-linear system, the two cases that follow the patterns (a) or (b) are actually 
trivial, in the sense that the field (j>\ can be seen as an artificial degree of freedom introduced in order to write the 
equation of motion of a (A^ — 1) fields as an A^-component coupled system. Case (c) is a bit more subtle. In fact, 
the Lagrangian of a irrotational barotropic fluid follows this pattern. We will work out this example explicitly in 
section 

Having completed the reduction, the new system of equations (jsj) can be analyzed as to its causal and geometric 
structures along the same lines as before. The new matrices f^'^ ab will determine the characteristics of the reduced 
system. Although the reduced system could now fail to be self-adjoint, the previous analysis can nevertheless still be 
applied — this is because (once spurious fields have been eliminated) the characteristics depend only on the leading 
symbol of the system of PDEs and are insensitive to lower order terms; cquivalently the Fresnel equation determined 
form the cikonal approximation is insensitive to any possible failure of self-adjointness. Once the behaviour of the 
fields (t>\ has been determined, the remaining spurious field (t>\ can (if desired) be obtained from by differentiation. 

F. Hidden approximate geometries 

Another situation of considerable interest occurs when, starting from a system of A^ coupled PDEs, one can neglect 
some of the degrees of freedom in some particular limited regime. In this case one can again perform a reduction 
process, but it is now important to remember that the reduced system is only an approximation to the exact behaviour 
of the system in the particular regimen. 

Let us take units such that time and length have the same dimension (some convenient reference velocity in the 
system is set to be 1) and the elements of Kab are non-dimensional. Then, the elements of must have the 
dimension of length, and those of the matrices f^'^AB the dimension of length squared. 

We can now assign to each dimensional coefficient in the set of second-order differential equations ([ill) a length 
scale. For instance, thinking of case (c), we might have a situation in which the length scales associated with the 
terms /^''n, T^^, f'^'^u, and f^" ai are all much smaller than those of the remaining terms in the system of equations. 
Let Ls be the largest length scale associated with the terms /^"n, f^^ib, and f^^ai, and let Li be the smallest 
length scale among all the other terms. Then, provided Lg Li, we can (in dealing with the physics at length scales 
larger than L^) safely approximate 

^"11,^,, r-'u, r\i^o (64) 

In these circumstances, for momenta low in comparison with 1/ Lg, (large wavelengths in comparison with L^) we will 
find that the fields 4)\ will behave as if they were coupled to an approximate effective metric determined by 

Tab = Tab - ^ {Kl ^b + Kl ^tb) (65) 

This is exactly what happens with the effective Lorentzian metrics arising in acoustic phenomena in dilute gas Bose- 
Einstein condensates (BECs) [isj. In this case Lg ~ h/ (to cq), (where h is Planck's constant, m is the mass of the 
atoms making up the gas, and cq is the velocity of sound in the condensate). This means that the geometric acoustics 
approximation in BECs is valid for wavelengths larger than h/{m cq) (an "acoustic Compton wavelength") at which 
one would probe the discrete nature of the condensed gas. 

The characteristic surfaces that one can construct based in these low-momentum metrics are only an approximation, 
and will only make sense when exploring the system with momenta lower that the scale \/Lg. By increasing the 
momentum we will find the true characteristic surfaces of the system. For example, in a BEC the true characteristics 
are those of a diffusion equation allowing infinite speed propagation. Nevertheless the effective metric found in low 
energies is a well defined Lorentzian metric with a maximum propagation velocity [ p^ . 

G. Fermionic fields 

The last possibility we shall comment on in this section is the extreme case in which all the matrices f^^AB are 
strictly zero. When this happens we cannot directly apply the previously derived analysis in terms of characteristic 
surfaces. What we should do instead is to either develop a modified discussion based on definition (|2^), or more 
prosaically, to first produce a system of second-order PDEs by iteration. Apply the first order operator 

Dab <^^^\ (r^B 9^0f + ^^(f:;^ <^f )) + Kab c^t (66) 
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twice. Then 

Dac Dcb c^i = 0. 
With this procedure we arrive at an equation of the form 

[rAB a.0f ) + f% a^0f + Kab <Pi = 0, 

with coefficients 



Tf^AB = Kac r^cs + Kbc r^cA + ^^,{^^'AC r^cB + T^Ac r^cs) + id^T^Ac) r^'cB, 
Kab = r'-Ac a,[(a^r^ci3) + Kcb] + [(d^r'^Ac) + Kac] [(^^r^cs) + Kcb]- 



(67) 
(68) 

(69) 

(70) 
(71) 



The system of equations ( pq ) now matches the pattern of our previous discussion. 

This situation is characteristic of the existence of fermionic fields. For example in the trivial case of the Dirac 
equation in flat space, we have T" ac = AC with AB now denoting spinorial indices. It is easy to see that 



AB = ^(t^ac l^CB + AC I^cb) = Sab r]^" , 



(72) 



so that the different spinorial fields all feel the same Minkowski metric ry'^". 



V. TWO-FIELD SYSTEMS 



If we are dealing with a simple two-field system, much of the algebra simplifies. This makes for a useful illustrative 
example. We start by writing Q{q,p) in terms of a completely symmetric four index tensor, (a quartic) 



Q{q,p) = Qiqr"" Pp P- = {F'niq) /''"22(g) - ^''12(9) M} Pp Po 



(73) 



In discussing the most general causal structure for the propagation of a two-field system it is this quartic QUiY'^p'^ 
that is the main geometrical object; instead of the metric f{q)^'^ — \J—g as in the single-field case ||l|. We 
can think of F^ ab as a 2 x 2 matrix in field space, f^" ab, whose components are point-dependent tensors (tensor 
densities, in fact): 



J AB 



F'll F\2 
^""12 F''22 



By redefining the linearized fields we can transform 

J AB 



J AB — l^A l^B J CD = AdAB J AB- 



(74) 



(75) 



In general, there is no combination of the three independent tensor components of f'^'^ ab, even with point-dependent 
coefficients, that vanishes. That is, typically 



Mab'\x) F\i{x) + Mab^^x) F''22{x) + 2Mab^\x) F\2{x) ^ 0. 



(76) 



In this algebraically most general case, we are forced to deal with pseudo-Finsler geometries. There are other cases 
in which the three tensor components of f'^'^ ab are not independent. Then, by using field redefinitions we can arrive 
at five different canonical cases: 

I: Two independent components: 

la: A diagonal component is set to zero — 



TAB 



fll fl2 
fl2 



(77) 
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In this case 



(78) 



lb: 



This case corresponds to quadratic reducibihty for the function Qiq^p) without diagonalization of the 
kinetic terms. There are two induced metrics, both based on = V~9 9^'^ ^ which identical to each 

other. However these induced metrics, while they successfully reproduce Q{q,p), and so reproduce the 
characteristic surfaces, are not enough to reproduce the symbol of the system of PDEs. From the induced 
metrics we can reconstruct Q{q,p) but not f^'^ABi so that we cannot arrange to diagonalize the system of 
PDEs in field space. 

The off-diagonal component is zero — 



Iab = 



hi 

/22 



In this case 



Q{q.p) = [r 11 P^.Pu) iF''22Pf.P.) 



(79) 



(80) 



It is this case that corresponds to a proper bi-mctric theory with diagonal kinetic terms and two distinct 
metrics. 

II: There is only one independent component: 

Ila; Only one diagonal component is different from zero — 



AB 



hi 





(81) 



In this case Q{q,p) = 0; this is a sign that one of the fields is unphysical, either because it is superfluous 
or because it is actually a gauge degree of freedom. 

lib: The two diagonal components are different from zero — and equal 



Jab = 



hi 






hi 



This is a true single-metric theory with diagonalizable kinetic energy terms and 

Q{q,p)^-{r\iP^.Puf ■ 
This is the case compatible with strict adherence to the Einstein Equivalence Principle. 
Ill: All the tensor components of Jab vanish: 



(82) 



(83) 



This probably means you are dealing with a Fermi field and somehow did not notice. See section (IVG) above 



In the three cases Ia^ h and lib the quartic ( |73|) factorizes into the product of two quadratics. The case /& corresponds 
to a proper bi-metric theory in which there are two fields each reacting to a different metric. In the cases la and //;, 
the two fields feel the same metric and so we would recover the usual behaviour of classical fields (such as Maxwell or 
Dirac fields) in general relativity. Of course in all these cases it is still necessary to analyze if the relevant metrics are 
Lorentzian. (And to verify that they satisfy suitable causality constraints.) The case Ila is more tricky: The quartic 
( [73| ) is degenerate and so we cannot strictly define any kind of characteristic surface. However, as we have previously 
argued, in some circumstances it is still possible to reduce the number of fields and so define characteristic surfaces 
based in the existence of a hidden metric. We shall now demonstrate this and related phenomenon explicitly with a 
pair of simple examples. 



VI. EXAMPLES 



A. The barotropic irrotational inviscid fluid 

We can write the Lagrangian for a barotropic inviscid irrotational fluid as |l^, 



^ = \p{^of + pdte + J^' dp' h{p'). 



(84) 
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Here, the two fields (j)"^ = (p, 9) are the fluid density p and the velocity potential 9 (the velocity field v can be written 
as V = 'V9 because the irrotational nature of the fluid considered). The function 



h{p) = h[pip)] 



dp' 



(85) 



lo Pip') 

is the enthalpy of the fluid. Now, by linearizing we obtain the relevant second-order Lagrangian for the perturbations 



-C^'^ - IpoC^Oif + piV9i ■ V9o + pidt9i + 

Z 2pQ dp 



pi- 



rn 



(Remember that the linear terms cancel by using the background equation of motion). From this we get the two 
linearized equations of motion 



(87) 



dt9i + vo • V^?i + ^ Pi = 0, 

Po J 

dtPi + V(poV6ii) + V(pivo) = 0, 

in which we have used the definitions vq = V^o and Cg = {dp/dp)\o- (Of course, these must and do agree with the 
equations obtained by the more traditional means of first finding the full equations of motion and then linearizing.) As 
we can see, equation ( ^7| ) has only the zero-order term for pi and in the discussion of hidden geometries the condition 
(c) above is fulfilled. In our previous notation (note that p and 9, which we boldface for clarity, are now field indices 
not spacetime indices) 



J f 



J pB = J 9p 

: 



0, 



60 







and 



while 



Kee = Kp0 = Kgp = 0, 



Po 



(89) 
(90) 



(91) 
(92) 
(93) 

(94) 
(95) 



Finally, by looking at case (c) as given in equation (|63|) we can calculate all the terms in equation (54). We have 
K = 0, T = and a single reduced metric of the form 



Jfiu ^ P0_ 
Co 



-1 



(96) 



which corresponds to the usual acoustic metric pg|, |2(]| . Once we have the solutions for 9i, the solutions for pi can be 
found by substituting in (^7|). 

On the other hand, if we tried to eliminate 9i to find an equation from which to derive the same set of solutions 
for pi, we would find a very involved equation. Formally one can write 



which implies 



Pi^T%ed^9i^V9i, 
9i^V-'pi. 



(97) 
(98) 
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Then we have 

5^(r"5.0i)-O, (99) 

and 

vOf^if'^d^v-^pi) = 0. (100) 

While at low momentum the differential equation satisfied by pi is rather complicated, we can see that in the large 
momentum limit (the eikonal approximation) the propagation of density waves sees the same causal structure as the 
^i-waves. That is, the causal structure for the propagation of both fields is described by the same metric determined 
byr^ = y^5^''. 

Of course with hindsight the reason this elimination procedure works is that it is possible to go all the way back to 
the Lagrangian in equation (^4|) and eliminate the density from the Lagrangian before linearizing. Nevertheless, this 
is a nice simple example of how to eliminate spurious fields from the linearized system. 



B. The BEC analog model 



Let us consider now a Bose-Einstein condensate. It can be described by the Gross-Pitaevskii equation, 

d 



dt 



ih-Mt,^)^ _ V^ + V;xt(x)+A |V'(t,x)r U(i,x), 



2 m 



where ^{t,K) represents the mean field wave function. By using the Madelung representation pi 



V'(t, x) = p{t, x) exp[—im6{t,x)/h], 
one can separate the (complex) Gross-Pitaevskii equation into the two (real) equations: 



dte + ^ivey + ^ + -p 

2 mm 



2m^ \y J 

dtp + v ■ {p ve) = 0. 



0. 



Here, the quantity 



Vq{p) 



2m \ ^ J 



(101) 

(102) 

(103) 
(104) 

(105) 



is what is commonly called the "quantum potential" [g2|, Expanding around a background solution of these 

equations one finds the following two equations for the perturbed quantities pi and 0i 

cl 



dtOi - vo • V9i - -^pi 
Pq 

dtpi + V{paV(f> 



TT^ D2P1 = 0. 
2m^ 



V(/9iVo) - 



(106) 



(107) 



We have used the definitions vo = V0o and c§ = (A/m)po- Additionally, D2 represents a relatively messy second-order 
differential operator obtained from linearizing the quantum potential. Explicitly: 



D2 pi 



1 -3/2 r« , +1/2ni , 1 -1/2 . , -1/2 ^, 

-2P0 [^(Po )\Pi+ ^(Po Pi) 



(108) 



These linearized equations are formally equivalent to the ones obtained for the barotropic irrotational fluid except for 
the presence of the quantum potential. This difference amounts to the existence of the following additional non-zero 
coefficients: 



pp 



2m^ 







2m2 



2m2 



0, Po'^'d^Po 

+1/2 



1/2 



-3/2 
-Po 



\i I -3/2^ +1/2 _ +1/2 

^(Po )+Po VPo -^Po 



(109) 

(110) 
(111) 
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Now, followi ng th e discussion in the "Hidden approximate geometries" subsection [V F , one can make the K coefficient 
in equation ( |l06| ) non-dimensional by multiplying the whole equation by Po/cq 



^Q. Then, the important point from our 
example is that it is easy to see that the new terms that occur in the BEC system (with respect to the previous 
fluid system) all come naturally multiplied by the the length scale — h/{m cq). This corresponds to an "acoustic 
Compton wavelength" defined by the propagation speed of phonons in the condensate. If one is probing the BEC 
system at l ar ge le ngth s cales (low momenta) in comparison with this Compton wavelength, one can neglect the new 
terms ( |109| ), (110), and (111). One can then perform the same trick as with the barotropic irrotational fluid to obtain 
an effective metric for pi. If we don't neglect the new terms coming from the quantum potential we can still arrive at 
an integral-differential equation for pi . This equation will encode the Bogoliubov dispersion relation ||l^, |l^ , showing 
us that the symbol of the Gross-Pitaevskii equation is that of a (complex) diffusion equation. As such it allows the 
propagation of signals at arbitrarily large speeds [p6|. 



VII. SUMMARY AND DISCUSSION 



Taking an arbitrary system of hyperbolic second-order PDEs (either generic or arising from a flrst-order Lagrangian) , 
the behaviour of the perturbations of the fundamental fields around any background configuration can be given a 
geometrical interpretation. For a single field there is always a nice and clean geometrical interpretation in terms of 
the d'Alembertian operator in an effective (typically curved) Lorentzian geometry Q. For several coupled fields, as 
discussed in this article, the situation is seen to be more complex: In simple cases all the fields see the same effective 
Lorentzian metric; with a little less luck the fields see different effective Lorentzian metrics, up to one effective metric 
per field; with extremely bad luck you will need to adopt some Finsler-like extension of the notion of Lorentzian 
geometry. Taking into account this formal analysis, in order to obtain an analog model of general relativity from a 
general system of PDEs it will be necessary to first require that the system fulfill some sort of "Einstein Equivalence 
Principle" (at least approximately). That is, at low energies we would desire an (approximate) factorization: 

F'^AB - hAB r - hAB 7=5 5^"- (112) 

Then, and only then, will all the low-energy fields see (approximately) the same effective Lorentzian metric, as is 
experimentally implied by the Eotvos experiment. At this level, we have not found it possible to obtain the "Einstein 
Equivalence Principle" from more fundamental principles. 

In searching for a geometrical interpretation coming from a set of PDEs we have analyzed the physical content of 
the matrix Jab- Sometimes, however, this matrix is singular and does not provide a geometrical interpretation for 
the behaviour of all of the fields. In these situations one can still find such geometrical interpretation by extracting 
additional information from the F^as and Kab coefficients. An extreme example is provided by a Fermi system. In 
this case, all the are zero but we can find a geometrical structure coming from the T^ab- In other situations, 
it can happen that the system has some "spurious" degrees of freedom making difficult the search for a geometrical 
interpretation in terms of the f^^AB alone. In many situations we can eliminate these spurious degrees of freedom and 
arrive to a reduced set of PDEs with a well defined geometrical interpretation. Additionally, there can be situations 
in which the elimination of spurious degrees of freedom can be justified in an approximate way. By this we mean that 
even though, strictly speaking, the Jab are non-singular and one should find a geometrical interpretation directly 
for the complete set of them, nevertheless for some particular regime (energy scale) the Jab can be considered to be 
approximately singular and an elimination of spurious degrees of freedom appropriate. The approximate geometric 
structure found by this procedure, though extremely useful, can be completely different from the exact geometric 
structure based on the exact characteristics. We have illustrated this point with the BEC analog model. The exact 
underlying geometric structure provided by the non-relativistic Gross-Pitaevskii equation is "parabolic" while the 
approximate effective geometric structure is Lorentzian or "hyperbolic" . 

In this paper we have discussed the different possibilities one can encounter in extracting a geometric structure 
from a system of PDEs. The particular systems that fulfill the Einstein Equivalence Principle can most easily be 
seen as analog models of General Relativity. Without additional structure they should not be thought of as models 
for general relativity, as the dynamics of the effective geometry can differ greatly from the proper general relativistic 
dynamics (the Einstein equations). On the other hand, in a previous work |jl[, and the in the context of a single field 
system, we showed that one-loop quantum effects could also provide, in some circumstances, a dynamics somewhat 
resembling that of general relativity. We leave as a future project to analyze in detail this quantum mechanism in the 
context of a multi-field system. 

In summary: The occurrence of effective metrics and effective geometries in low-energy linearized approximations 
to wide classes of dynamical systems is striking — the near ubiquitous occurrence of this effect is particularly useful 
when considering analog models of general relativity, can even be interpreted as being strongly suggestive that some 
form of "induced gravity" may underly true physical gravity. To turn "induced gravity" into a serious contender 
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several additional conditions must first be met: On the one hand we will need to be dealing with multiple background 
fields (at least six) in order to have an algebraically general spacetime metric [Q , on the other hand all the linearized 
background fields must see the "same" Lorentzian metric (either exactly the same metric if you believe strongly in 
the Einstein Equivalence Principle, or at the very least approximately the same metric in order to be compatible with 
experimental constraints from Eotvos-type experiments). Only once these two basic conditions have been met will 
it be sensible to develop a Sakharov-like "induced gravity" ||2J|. One-loop effects will then generate a term in the 
effective action that is proportional to the Einstein-Hilbert action jl] . The result would then be an "effective" theory 
of gravity in the sense of DonoghuejE5| , but would still suffer from the potential defects common to all embedding 
models of general relativity ||l|, |27|. 

In this regard the results we report in this paper are mixed: The good news is that the occurrence of Lorentzian 
metrics seems quite common, the bad news is that multiple Lorentzian metrics (and worse) seems generic. We have 
been able to say quite a bit about the general way in which one might probe the causal structure of these theories, and 
while we have been able to put much of the discussion of multi-refringence in general framework. We are continuing 
to work on seeking field theories that are more general than simple "minimum substitution" theories, but are still 
restricted enough to produce a well-controlled metric structure. 
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